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ABSTRACT 


Advances in nuclear and diesel-electric submarine 
technology have reduced the effectiveness of passive means of 
detection. The United States is faced with a multipolar 
mareat im part due to the proliferation to Third World nations 
of advanced diesel-electric submarines. The use of active 
sonar must be explored to gain back the detection advantage 
the United States submarine force has enjoyed in the past. 
The use of bistatic sonar reduces the counter-detection threat 


mesulting from active sonar. 
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THESIS DISCLAIMER 


The reader 1S cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification 1s at the £2sieert veneer 
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I. INTRODUCTION 


A. PURPOSE 

The purpose of this thesis is to study an application of 
active bistatic sonar to anti-submarine warfare (ASW). In an 
active monostatic sonar system the source platform acts both 
as a source and receiver. A bistatic sonar system 
incorporates a separate source and receiver platform. 
Specifically, this research will concentrate on the use of an 
Mi@tave bistatic sonar system to protect a channel from a 
transiting enemy submarine. A mathematical model can be 
developed to calculate signal excess as the enemy submarine 
transits down the channel for a given location of the source 
and recelver. This signal excess can be represented as a 
detection probability, which will be used as the payoff ina 


Mmere-DerSon zero-Sum game. 


B. BACKGROUND 

Recent political changes in Russia and other former 
communist countries have resulted in a significant reduction 
of the perceived Eastern Bloc threat to the United States. 
The immediate cost of the Cold War victory is decreased world 
Stability. Diminished Russian global influence has caused the 
Smstdbility throughout the Third World [Ref. l1:p. 146]. 


Reduced world stability has forced the United States to face 


a new Mlttpolar Sehmeaee Consequently, U.S. milit@aa 
strategists now place a much higher priority on reg vomem 
Cont lice plammiunge For the ASW mission, the United States 
must now also contend with diesel-electric submarines 
operating on their batteries, quieter nuclear submarines, and 
a shift in policy from preparing for global wartarpeumee 
littoral warfare. The Russian navy, especially their 
submarine force, remains a viable fighting force. 

The ASW advantage enjoyed by the United States and its 
allies has steadily decreased. Over the last ten years 
Russian submarines have undergone significant reductions of 
their radiated noise signatures [Ref. 1:p. 146]. During ene 
same period improvements of passive detection by the United 
States advanced at a slower rate. The proliferation of 
advanced diesel-electric submarines to Third World nations 
adds a new dimension to the ASW problem for the United States. 
The new diesel-electric submarines are capable of extended 
patrols without having to snorkel as frequently as in the 
past. New battery technology provides a cheap alternative to 
nuclear submarines for Third World countries. To maintain its 
advantage over other countries in ASW, the United States must 
revolutionize detection technology. 

1. Historical Perspective 

The early 1960s saw the beginning of the tactical and 
strategic development of nuclear submarines. The United 


States development was significantly different from the Soviet 


Union’s. The United States invested in maintaining vital sea 
lanes of communications and assuring freedom of navigation in 
forward areas. To support this global ASW mission, the United 
States concentrated on the development of highly capable 
passive sonar systems, gaining a significant ASW advantage. 
For three decades after World War II the United States 
submarine force enjoyed a tactical and strategic advantage 
over the Soviet Union’s submarine and ASW forces. 

Since the late 1970s the United States ability to 
detect enemy submarines by passive means has significantly 
decreased. During the period 1975 to 1988, the radiated noise 
Signal levels of Soviet submarines dropped by 30 dB or by a 
factor of 1000 [Ref. l:p. 146]. The increased use of diesel- 
electric submarines by Third World countries compounds the 
passive detection problem. 

Active sonar does not rely on target generated noise 
as with passive detection. The greatest concern for active 
sonar 1S counter-detection. ft =e very  drificult for a 
Submarine to remain undetected when using active sonar. 
However, accurate bearing and range information can be 
gathered by using active sonar. With passive systems only 
bearings can be found, thus requiring target motion analysis 
to find the range to the target. 

2. Bistatic versus Monostatic Active Sonar 
The most significant problem of active sonar is the 


loss of stealth of the searching platform. An active sonar 


can typically be detected by another platform at twice the 
detection range of the searching platform. The sound mug 
travel both directions for the searcher but only one direquauem 
for the potential target platform. 

With the use of a bistatic active sonar, where the 
transmitter (source) and receiver are two different platforms, 
the searching ship can remain undetected but also take 
advantage of the active sonar. It 1s, however, important to 
ensure the receiver 1S not counter-detected by the target as 
active sonar does not discriminate between searcher and 


target. 


C. DESCRIPTION OF BISTATIC SEARCH SCENARIO 

A source and receiver platform are searching out a 
rectangular barrier area attempting to detect and locate a 
transiting enemy submarine. The target submarine’s 
Orientation is known by the receiver and source platform. The 
receiver is placed within the barrier for two reasons. First, 
to constrain the problem to help make it solvable, and second, 
to keep the receiver within range to prosecute the target if 
necessary. The source platform can be located inside or 
outside the barrier area. The target enters the barrier along 
one side and transits through. Figure 1 shows the target 


etry jAnto the barriers 


Barrier Area 


Target 





Figure 1. Target entry point into the barrier 


The detection probability can be calculated, for a fixed 
source location, for each given case of target entry and 
receiver location. These probabilities represent a payoff 
matrix. On one axis of the matrix are the receiver locations 
and the other axis are the target entry points. To determine 
the result for a given combination of receiver location and 
target entry point, the detection probability for the receiver 
1s the intersection point within the payoff matrix. The 
receiver does not know where the target will enter the 
barrier. The target also does not know where the receiver is 
located within the barrier. The receiver desires to maximize 
this value of detection, where the target wants to minimize 
the value. This method of one side maximizing and the other 
Side minimizing 1S a two-person zero-sum game [Ref. 2:pp. 10- 


ae. 


There are two possible outputs to a two-person zero-sum 
game, a pure strategy and a mixed strategy. The pure strategy 
means that the element within the matrix represents both the 
largest column and the smallest row value [Ref. 2:p. 11] or in 
other words, the matrix contains a saddle point. The mixed 
Strategy occurs when no saddle point exists, thus there is a 
probability distribution of pure strategies [Ref. 2:pp yee 
dele 

In the evaluation of the payoff matrix formed by the 
model, the underlying desire 1s for a single pure strategy for 
any given case of a fixed source location. If there dem 
saddle point, the model results provide the receiver an 
optimal pure strategy to detect a transiting submarine through 
the barrier. If however, the results yield no saddle point, 
a mixed strategy 1S required and thus the receiver and the 
source must make tactical decisions regarding the choice Gf 


strategy. 


Il. MATHEMATICAL MODEL 


A. MODEL ASSUMPTIONS 

The mathematical model incorporates the bistatic active 
sonar equation to determine the optimal strategy for a 
receiver and a source to detect a transiting enemy submarine. 


The bistatic active sonar equation is 


Si =i ieee ea NE ne TS HAG. (1) 

Where SE is the mean signal excess, SL is the source level of 
the active transmission from the source platform, TL,, 1s the 
transmission loss from the source to the target, TL, 1s the 
transmission loss from the target to the receiver, NL is the 
ambient noise level of the surrounding sea, DT is the 
detection threshold of the receiver sonar system, TS is the 
bistatic target strength, and AG is the array gain of the 
receiver sonar. All parameters are measured in decibels (dB) 
and referenced to 1 micropascal. 

Signal excess (SE) is assumed to be a normal random 
variable with mean given in equation 1 and standard deviation 
Oo. Detection occurs whenever SE>0. ihe probabilicy “or 


meeection 1s 





P(SE>0) =© (52), (2) 


where ® is the cumulative normal distribution and o is the 


standard deviation [Ref. 3:p. 3-2]. Typical values Eor )G sane 
between 3 and 9 dB [Ref. 3:p. 3-2]. The value of Oo for Ehe 
model 1s set at §~S “de: For computational purposes, a 
polynomial approximation to @® from the Handbook of 
Mathematical Functions [Ref. 4:p. 932) waseuces 

The fictional active low frequency source is assumed to 
have a source level (SL) of 220 dB re 1 micropascal. This 
source level provides a theoretical detection range in excess 
of 300 nautical miles (nm) [Ret. Peo eee 

Transmission loss (TL) results from attenuation and 
spreading. The model assumes the spreading loss is due to 
spherical spreading with no transition to cylindrical 
spreading. This assumption yields reasonable fits to actual 
measured data under various conditions, Urick [Ref. 5:p-) aia 


Transmission loss can be expressed as 


TL=20 logR+aRx10°7, (3) 
from Urick [Ref. 5:pp. 110-111], where R is the range in yards 
from either the source to the target or the target to the 
receiver and @ is the absorption coefficient im eewie 
ka lbovardy (kya For the model @ = 0.0041 dB/kyd andi 
CaleuLabed. from 


se Oe ee ee 


—— es 75 x1 ee (4) 
L+£* <4 .GnOO a 


a 


where f£ 1s the search frequency expressed in kHz [Ref. 5:p. 


108). For the model, f is assumed to beaut eie 


The background noise 1S assumed to be ambient noise 
limited, vice reverberation limited. The average noise level 
(NL) 1s 63 GB re 1 micropascal for the search frequency of 100 
Hz found on the Wentz-Knudsen curves [Ref. 5:p. 210]. Figure 


2 shows the Wentz-Knudsen curves. 
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Figure 2. Average deep-water ambient-nolise spectra. 


The detection threshold (DT) for the receiver sonar is 12 
adB re 1 micropascal. Detection threshold can be calculated 
from one of two different cases [Ref 5:p. 394]. Case I 
assumes a known signal in Gaussian noise and a fixed point 
target at a known range. Case II assumes an unknown signal 
and a target of unknown range and range rate. Case II yields 
a higher value and thus 1S a more conservative estimate of 


detection threshold. Case II detection threshold is 


se 


(ayes jreye; == (5) 


where d, the detection index, is determined by the specified 


probability of detection and probability of fEalse Walon 
Selecting values for probability of detection (Pd) equal to 
90% and probability of false alarm (Pfa) equal to 0.01%, 
yields d=25. Wis the frequency bandwidth (100 Hz) and t is 
the pulse duration (10 seconds). [Ret >) seo 

The bistatic target strength (TS) is a function of the 
minimum (TS,,,) and maximum monostatic target strength (2ouam 
and the bistatic angle (0,) formed by the relative angles, (@,) 
and (@,), between the source and the target, and the receiver 
and the target, respectfully. Figure 3 contains a Viste 


representation of the bistatic angle. 


Transmitter 


6 ,= 0.5(6R + 87) 


Receiver 





Figure 3. Bistatic angle for target strength calculatvaas 
The equation for target strength is given by 


TS= (TS..,-7S 


m 


in) > DS aS ene (6) 
where, for a typical submarine, TS,,,=25 dB and TS:,,=12 dBure 


l-micropascal. [ReEt. 6:ponZeole 


ie 


Taeweimal term for the active bistatic sonar equation is 
the array gain (AG). Array gain is a measure of the sonar’s 
ability to increase the ratio of the submarine signature to 
the background noise. Array gain of 15 dB re 1 micropascal is 
consistent with the conclusions of Stefanick’s work [Ref. 7:p. 


me |. 


B. COMPUTER MODEL 

The code for the computer model is written in Turbo Pascal 
6.0 for DOS based computers. The computer program calculates 
Mme detection probability of the target as it transits through 
the barrier for a fixed entry point Fore (ane target anda fixed 
receiver and source position. 

1. Inputs to Model 

There are two input files for the main program. The 

mimest file contains specific information about the barrier, 
receiver, target, and source. The dimensions (width and 
depth) of the barrier are provided by the input file. Inputs 
for the receiver consist of the step size for receiver 
location in both the X and Y axis, and the detection threshold 
and array gain of the receiver sonar system. For the target, 
the step size for the target’s barrier entry, speed of the 
target during transit and its minimum and maximum monostatic 
target strength are provided. The source inputs include 
active sonar pulse frequency, pulse length, ping interval, 


source level of active sonar and ambient background noise 


ey 


level. The second input file contains all source locations to 
evaluated. 
2. Main Program 

The program calculates for each given receiver 
poSition and target entry point, with a fixed source positvom 
the detection probability for target transit through 
barrier. The target entry points are based on the user input 
of the width and depth of the barrier, and the step size for 
target entry. For example, if the barrier is 200 kiloyards 
(kyds) in width, 40 kyds in depth and the target entry step 
size 1s 10 kyds then the target would enter the barrier at X 
position equal to (0,10,20,30,...,200) kyds with Y posi 
equal 40 kyds. Based on the same width and depth of the 
barrier as above and with user input of step size of 5 kyds 
for both the X and Y position for the receiver, the receiver 
would be located at all possible combinations of X positions 
of (0,5,10,...,200) kyds and Y positions of (075300; 
kyds. 

To determine the overall detection probability for the 
target transit, the signal excess is calculated using the 
active bistatic sonar equation for the receiver at fixed 
points for the target on the transit. The Signal excess can 
be represented as a detection probability as prevvousday 
discussed in the Model Assumptions section. The probability 
of no detection is found by taking one minus the probability 


of detection. By taking one minus the product of all the 


nz 


probabilities of no detection the overall detection 


probability is determined as shown by 


P,=1-II(1-p,), (7) 
where P, is the overall probability of detection of the 
Meameit and p, is the probability of detection of the target 
Bemene 2° ping. Each P, is assumed to be independent as the 
time between pings is typically greater than 5 minutes. The 
fixed points are based on the speed of the target, ping 
interval of the source active sonar and the length of the 
barrier. Fach ping of the active sonar represents a fixed 
Beant for the target. Between pings the target advances 
through the barrier based on its speed. 

Each possible entry point for the target is evaluated 


for detection probability (A,,,) for a fixed receiver and 


XYZ 
Mamas OSitIion. A, 1S the probability of detection for a 
given (x,y) coordinate position of the receiver and entry 
point Zz for the target. This procedure is repeated for all 
possible receiver position combinations. The results for each 
possible combination of target entry and receiver position is 
written to an output file. The program is rerun for each 
possible source location. The code is shown in Appendix A. 
3. GAMS Interface 
The output from the program forms a probability matrix 


with one axis receiver position and the other target entry 


moint. This probability matrix is for a given source 


ile 


DOSITE Lom. Fach different source position forms anoemew 
separate probability matrix. Using Game Theory, the 
probability matrix can be evaluated to find an Opti 
strategy for both the receiver and the target. To £Lindiieme 
optimal strategy numerous simultaneous equations must be 
solved. This problem can be formulated into a linear program 


as shown below: 


maximize V 
subject to } Yo A, P,,-V20; z€q, 
xemyen (8) 
P,, =1, 
xEmyen 
and Pia je, Eee 


where m and mn represents the number of different X and Y 
positions for the receiver respectfully, and ¢g is the number 
of different entry points into the barrier for the target. V 
represents the overall probability of detection of the target 
for the optimal strategy (value of the game) [Ref. 8:pp. 3l- 
32). A,, 1s defined in the previous section. The z value for 
the target represents the X coordinate position as the Y 
coordinate position is at the top of the barrvem PS 
represents the percentage the receiver should spend in a given 
(x,y) coordinate position. The typical problem size, used 
specifically for this model, had 22 constraining equations and 
369 variables. 

The General Algebraic Modeling System (GAMS) MINOS5 


linear program solver, version 2.25, was used to solve the 


linear program. Appendix B contains the model file for the 
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solver. DOS batch files were used to integrate the Pascal 
code and the GAMS solver. 
4. Counter-detection Determination 
Three methods were considered for calculating the 

counter-detection threat. These methods correspond to three 
assumptions about what happens after a counter-detection: 

me nothing, 

2) target evades, and 

3) target attacks receiver. 
The target’s detection threshold is higher than the receiver’s 
as the target does not have the same sonar system or the 
characteristics of the active pulse from the source platform 
as the receiver. For simplicity, the target’s sonar system is 
assumed conservatively identical to the receiver’s' for 
counter-detection calculations. The following notation is 
used ee) help define the various counter-detection 
calculations: 

I = ping on which receiver detects target for 1° time, 

J = ping on which target detects receiver for 1°° time, and 

Ior J=o, if no detection during a given time period. 


With this notation the P, 1n equation 7 is the equivalent of 


P(I<m), (9) 


where mis the time required for the target to transit through 


Pie barrier. 


ie: 


The first method determines counter-detection three 
by calculating a counter-detection probability and comparing 
the result to the detection probability. Counter—-detecuaas 
probability is based on the optimal strategies determined by 
the linear program solver. The program, based on the optimal 
receiver location and target entry point, determines the 
counter-detection probability of the receiver. The counter- 
detection probability is multiplied by the probability the 
receiver will be located in that position and the probability 
the target will enter the barrier at that point. All possible 
combinations of optimal positions for both the target Wane 


receiver are summed up as shown by 


P.g=& [1 -II (1 -P.a) ] P,,P,. (10) 
where P., represents the overall counter-detection probability 
of the receiver for a given source location. P., 1S jae 
counter-detection probability of the receiver on the 2°79 pumas 
P,, 1s the mixed strategy probability for the receiver amemee 
1s the mixed strategy probability for the target. The mixed 
strategy probability for the receiver comes from the solution 
6f the linear program in equation sce The mixed strategy 
probability for the target is the dual solution of the linear 


program in equation 8. 


lis 


The second method determines the counter-detection 
threat calculating the probability of the receiver detecting 
the target before the target detects the receiver and before 


time elapses, shown by 


P(I<J and I<). (11) 


This takes place of the detection probability of the target, 


P(I<~) =P(receiver detects the target), (12) 


in the payoff matrix. 
The third method finds the probability of the receiver 
detecting the target first and subtracts the probability of 


the target detecting the receiver first, 


Pilko) =PCrc Er, (13) 


and the result is stored in the payoff matrix. This method 
will show if the target will have a significant chance of 
detecting the receiver before the receiver detects the target 
and thus be able to attack first. 

The payoff matrix for the last two methods, when 
evaluated using the linear program solver, yielded a strategy 
for the receiver very similar to the original method one 
payoff matrix. For model run-time considerations, only one 
method of counter-detection threat assessment is used. The 
first method is incorporated into the model as the strategy 
results are similar to the other methods and a specific value 


mem counter-detection is found. 


de 


5s Output 

The output from the GAMS linear program Solver 
written to an ASCII file and a Pascal program converts the 
information to the final output file. Appendix B contains eae 
code for the conversion and output to the final output fae 
This allows the output file to be appended as multiple 
independent runs must be done to include data from all 
possible source locations. The final output file contains the 
optimal strategy for the receiver and the target, the overall 
detection probability (or value of the game), and a counter- 
detection probability of the receiver for each given source 


LOocat Lom. 
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III. RESULTS 


A. INTERPRETATIONS OF RESULTS 

The model uses the sonar equation assumptions discussed 
previously. The values for the barrier dimensions, receiver 
and source locations, and target entry points are also 
required for the model. The barrier dimensions are set at 200 
femlovyards (kyds) in width and 40 kyds in depth. The target 
enters the top of the barrier at 10 kyds intervals and 
proceeds through. The possible receiver locations, within the 
barrier, are at intervals of 5 kyds in both the X and Y axis. 
Figure 4 graphically shows the barrier, target entry points 
and possible receiver locations. 
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Figure 4. Target entry points and receiver locations. 
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Figure 5 shows the barrier imposed on all 4,469 possible 


SOuLece Vike eat ite mom 


Possible Source Locations 


Source EE 
EE Borrier 


J 


300 


oe the ori 0.0) ts locat the 
X Position in kiloyards eS Coe Nae 


“ 
Uv 
i= 
o 
> 
= 
“x 
= 
c 
= 
— 
a 
Oo 
oe 
> 





Figure 5. Possible source locations. 


Note, in Figure 5 what appears to be straight vertical lines 
are actually many small circles, representing source 
locations, very close together. 

At the completion of a model run for a given source 
location, an optimal strategy and the value of the game 
(overall detection probability) for the receiver is found. 
The dual solution of the results provides the optimal barrier 
entry point(s) for the target. In virtually every case the 
optimum strategy 1s mixed for both the receiver and the target 


platform. For example, with source located at the center of 


20 


mee barrier (100,20), the optimal pure strategies for the 


meceiver were as follows: 


xe Y ) Probability 


(45;,.00') 0234 
Go, 0) O07 
(Sor 20) Gs bhesS, 
0) 0.068 
(145,20) 0.224 
(Ss, ay Oe, LOG 


and for the target the optimal pure strategies for entry into 


mae barrier were: 


(Z) Probability 


(0) 0.437 
(30) oe 0lats: 
(100) OU So 
Clie) OG bs 
(200) Oman. 


In almost all other cases the optimum pure strategies were 
more numerous than for the example above. 

Figure 6 shows the overall detection probability (value of 
the game) for various source positions. 
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Figure 6. Value of the Game for various source positions. 
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Each of the 4,469 points in Figure 5 required separate solving 
on a computer linear program solver. The highest overall 
detection probability in Figure 6 occurs in two places asea 
source located just off the center of the barrier at xX 
position 100 kyds and Y posittt96m 20 en ee eee At tha 
location the receiver has a overall detection probability of 
G5 7% This high value for detection demonstrates Vere 
usefulness of bistatic sonar. The detection probabilities are 
symmetric about the center of the barrier on both the X and Y 
axis. Values of detection probability around the center are 
fairly close to the highest value. For this reason the 
surface in Figure 6 has a well-rounded top. The receiver is 
located at least 20 kyds away from the source, based on the 
optimum mixed strategy. In most cases the receiver is located 
in excess of 60 kyds from the source. 

The most significant problem of placing the source near 
the center of the barrier is the counter-detection threat to 
the receiver. Figure 7 shows the counter-detection 


probability for all possible source locations. 
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Figure 7. Counter-detection threat. 


The apparent noise of the counter-detection threat plot 
results from the calculation method used. Counter-detection 
probability is based on the optimal solution for the receiver 
to detect the target. The receiver avoiding detection by the 
target 1S not optimized. 

Counter-detection probability is reduced significantly by 
moving the source either up or down or left and right of the 
center of the barrier. With the source located off the center 
of the barrier, the detection probability is reduced only 
slightly. Figure 8 shows the contour graph of the counter- 
detection surface and better illustrates the sharp decline of 


@eumeer-detection probability. 
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Figure 8. Contour graph for the counter-detection threat. 


At the conclusion of a model run data provided includes 
overall detection probability, mixed strategy for the 
receiver, mixed strategy for the target and counter-detection 
threat to the receiver for a given source location. Multiple 
runs of the model with varied source position allows for the 
determination of the best strategy for both the receiver and 
the source. These results contribute a tool for decision 
making regarding the best receiver and source placement to 


detect an enemy submarine. 
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B. UTILITY OF APPLICATION 

The model provides a decision making tool regarding 
submarine detection. Many different possible source and 
receiver position scenarios can be evaluated with the output 
data. Emphasis on detection of the target or counter- 
detection of the receiver can easily be incorporated into the 
decision making process. 

The model provides a theoretical prediction of detection 
probability of an enemy submarine transiting through a 
barrier. The model results allow a user to evaluate the best 
strategy given any number of constraints. These constraints 
include source and receiver locations, barrier size, and 
farget entry point. Limitations on these factors could be a 
result of geographical considerations or some other physical 
constraints. 

Importance can be placed on either detection or the 
counter-detection probabilities by review of the output from 
the model. If the user 1s not concerned with being counter- 
detected more emphasis can be placed on detection probability. 
If however, the enemy submarine 1S very capable then the 
counter-detection threat can be considered with greater 
importance. These tools provided by the model allow the user 
to make logical decisions based on the current situation. 


Counter-detection cannot be minimized, however. 
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IV. PROGRAM 


A. SYSTEM REQUIREMENTS TO RUN MODEL 

The code portion of the model, written in Pascal computer 
language, can be used on any IBM compatible personal computer. 
The GAMS program requires at least a 386SXK IBM compatible 
personal computer to execute the linear program solver for the 
model. A math coprocessor 1S not required, however it speeds 
up both the Pascal code execution and the GAMS linear program 
Solver. Only 640 kilobytes of random access memory (RAM) 1s 
needed for execution of the model, but more RAM speeds up the 
GAMS linear program solver and allows for a virtual drive to 
be created. The virtual drive is used to speed up the storing 


and retrieving of temporary data during model runs. 


B. LIMITATIONS OF THE CODE 

The most significant problem with the model 1s execution 
Eime FOr {the -pregitam. The model ran on an IBM compatible 
486DX-25 MHz machine with 8 megabytes of RAM. Each run for a 
given source location took two hours to complete. When the 
model was run on a 386SXK-20 MHz machine with math-coprocessor 
and 5 megabytes of RAM, execution time was six hours for each 
run. The model run time factor limits the number of source 
locations to be evaluated. The model runs takes advantage of 


symmetry about the center point of the barrier as discussed in 
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the Results section. These long run-times are in spite of 
considerable effort spent in writing code efficiently. For 
example, the equation to determine target strength contains a 
Sine squared (sin‘*(8,)) term. With the use of some 


Mergonometric functions, sin*(®,) could be represented by 


il i 2 
Ropkrp Worlop + = [Ror - (Ror ~ Rp) Alie (14) 
[Ref. 9], where R,, 1S the distance between the source and 
earget, Rr 1S the distance between the target and the 


receiver, and R,, 1s the distance between the source and the 
receiver. The other two terms (h,, and h,,) are the horizontal 
distances between the source and the target, and the target 
and the receiver respectfully. The horizontal distances were 
the only new calculation required, as the distances between 
platforms are already calculated (each requires a square root) 


Mmemecne transmission loss terms. 
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V. CONCLUSIONS 


A. SUMMARY 

The use of active bistatic sonar provides a new method of 
detecting very quiet nuclear and diesel-electric submarines. 
Bistatic sonar takes advantage of the bearing and range 
information from active sonar and the stealth of passive sonar 
by using a separate source and receiver. The development of 
a model incorporating bistatic sonar allows for the 
determination of the optimal strategies for both the receiver 
ang “Source platloums. 

Multiple model runs are required to find the opal 
Strategy for the receiver and the source. In each case 
evaluated, for a given source location, the program yielded a 
mixed strategy for the recelver. Comparing the overall 
detection probability (value of the game) for each model run, 
for a given source location, provides a method for optimum 
source placement. There are multiple source locations 
yielding similar detection probabilities. 

The mixed strategy for the receiver, with a given source 
location, places the receiver at least 20 kiloyards (kyds) 
from the source and in most cases outside 60 kyds. Lf the 
optimum receiver strategy placed it very close to the source, 
the use of bistatic sonar would be the same as a monostatic 


SOltal. 


ae 


The counter-detection threat to the receiver is calculated 
based on the optimal strategy of detection for the receiver. 
The highest threat occurs within the middle of the barrier. 
By placing the source away from the center of the barrier 
counter-detection threat drops significantly while the 
detection probability remains fairly constant. 

All the previous results provide tools for the employment 
of a receiver and a source platform to best detect a 
transiting enemy submarine. With the counter-detection threat 
and mixed strategy results for the receiver, many different 
combinations of receiver and source locations can be 


@enmsidered. 


B. POSSIBLE FOLLOW-ON RESEARCH 

Two possible avenues for further research include: use of 
aera l target strength, sonar system parameters and 
environmental data; and simulation of mixed strategies for the 
receiver, target and source. The use of real data for target 
strength, sonar system parameters and environmental data, 
would yield more accurate results. It also allows for 
evaluation of an actual search scenario with a known searcher 
and target platform. This approach would make the model more 
of a tactical decision aid. 

Developing a simulation program would help confirm the 
results yielded from the model. ier would Galllow forthe 


comparison of the model to data generated by repeated 


BE: 


Simulations. AT&T Bell laboratories have 





Simulation program, Submarine Surveillance Simul 


that could be used in this capacity. 
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APPENDIX A 


A. MAIN PROGRAM 


Program Thesis; 

{ 

LT Aasgeir Gangsaas USN 

Naval Postgraduate School 
midesis work on Bistatic Sonar 
Turbo Pascal 6.0 


The main program reads in the user inputs about the 
receiver, target and the source. Also the user inputs about 
the environment. With this information the program 
calculates the signal excess as the target transits through 
the barrier. This signal excess 1S converted into a 
probability. The program contains three subroutines, 
ProperSourceData, ReadIn and DetermineProb that are 
explained below. } 


Uses DOS, Normal, Signal4, Gamsform; 


Sonsct Sigma = 8; 


Var Sx, Sy, WBarrier, DBarrier ; Real; 
Dz, Dx, Dy, Freq, PingIn, Num : Integer; 
PulseL, TSpeed, Rx, Ry, Tx, Ty : Real; 
iomax, TSmin, NL, SL, AG, DT >: Real; 
Done : Boolean; 


Ae ORR ae ERE REAR EEK KERR ERE KK E KK ERE | 


Procedure ProperSourceData; 


{ This procedure reads in the source location data for 
multiple source locations. It ensures the proper source 
location sequence for the multiple model runs.} 


far lt, J, N : Integer; 
hm ¥ >: Real; 
mono, Datalnscource : Text; 


BEGIN 
Assign(Temp, ‘D:\TP\WORK\Num.txt’) ; 
Assign (DataInSource, '’D:\TP\WORK\THESSRCE.IN’) ; 
Reset (DataInSource) ; 
Reset (Temp) ; 


ou 


Readln(Temp,J) ; 

Readln(DataInSource) ; 

FOR LL. . =) Sper ube 
Readln(DataInSource,N,xX,Y); 

IF not eof (DataInSource) THEN 


Done := False 
ELOE 

Done ~s= Truc 
Se Sa 
She — ae ae 
Nui SS aN 


Rewrite (Temp) ; 
Writeln(Temp, (Num + 1)); 
Close(Temp) ; 

Close (DataInSource) ; 


END; 


[EERE KEKKEEKKEKKEKKKK KKK KEK EKER EE EERE ES So ee 


Procedure ReadIin; 


This procedure reads in from a file (Thes.in) 


information about the receiver and the target needed for the 
program. } 


VaroDatalni=) Trexte 


BEGIN 


Assign (Datalin, 'D:\TP\WORK\THES.IN’); 
Reset (DatalIn) ; 

PropersourceData; 

Readln(DatalIn); 

Readln(DatalIn, PulseL) ; 
Readln(DatalIn) ; 
Readln(DataiIn, Freq); 
Readln(DatalIn); 
Readln(DatalIn, PingIn); 
Readln(DatalIn) ; 
Readln(DatalIn,WBarrier) ; 
Readin(DatalIn)  ; 
Readln(DataiIn, DBarrier) ; 
Readln(DatalIn) ; 
Readln(Datalm nz): 
Readln(DatalIn) ; 
Readln(DataIn, TSpeed) ; 
Readln(DatalIn) ; 

Readln (Datalin, Dx) ; 
Readln(DatalIn) ; 
Readln(DatalIn, Dy) ; 
Readln(Dataln) ; 
Readln(DataiIn, TSmax) ; 


Be 


Readln(DatalIn); 

Read lmDatalLn, Pomin) ; 
Readln(DatalIn) ; 
Readln(DataIn,NL); 
Readln(DatalIn) ; 
Readin(Datain, SL) ; 
Readln(DatalIn) ; 
Readln(DataiIn,AG) ; 
Readln(DatalIn); 
meadln(Datalin,DT) ; 

END; { Procedure ReadiIn} 


LGR Reece ears enh eee Ke ACR RK RR RK RR RK KK | 
Procedure DetermineProb; 
{ This procedure calculates the detection probability 


meee FOr the receiver £or each receiver location and 
target entry point.} 


wee VOtalDz, TotalDx, TotalDy : Integer; 
mbeece . €o> TotalTime, Ti T2 =< Integer: 
MmeMnp o5, Freadsq, a, tq : Real; 
Wawaout = Text; 
BEGIN 
ReadIin; 


Sesiogn(DataOut, ‘F:\THES.OUT’ ); 
Rewrite(DataOut) ; 


ieee = (). 0; 

ite: = DBarrier; 

oan — |), O- 

ame = 0.0; 

MetalTaime := Trunc( (DBarrier/ (TSpeed*2))*60) ; 

on := 0; 

eros = 0; 

a 2 = 0: 

FreqsSq := SQR(Freq/1000)  ; 

a:= (0.1*FreqSq/ (1+FreqSq) )+(40*FreqSq/ (4100+FreqSq) ) 


+((2.75E-04) *FreqSq) ; 
WHILE (C3 <= WBarrier) DO 


BEGIN 
Vite (G2 —=|} DBarrrver) DO 
BEGIN 
WHILE (Cl <= WBarrier) DO 
BEGIN 

ie) 

Whee Wile = TotralTime) DO 
BEGIN 


Sime solic xcesa (Sx, SV, Tx, Ty, Rx, Ry, 
Euubscin, a, TSmax, TSman, 
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NL, SL, AG, Die 
te Tl =O) eee 


Temp := 1 - NormalValue(SE/Sigma) 
ELSE 
Temp := Temp* (1-NormalValue(SE/Sigma) ); 
Ty :-= Ty - (TSpeed/(3* 20) Fimgie 
Tl. 2=. TLS Pingine 
EN: 
Ty #= DEarrrer- 
TZ 
Write (Dataout, (1 - Temp) eee 
Cl 73 =26 wb 
Te v= Tx +. bz; 


END: 
Writeln(DataOut) ; 
C2* = 36 ee 


Jip ea S10. 
Rye s= aol 
Cie -10- 
END 
C37 -=- C3. be: 
So. Cua — ee OS ae 
CZ a=) 20% 
Rye = eo, 
END: 
Close(DataOut) ; 
END; 
fA A ee ee ee Maan Program KERR EE EE AX * te 
BEGIN 
DetermineProb; 


FormatG(Trunc (WBarrier) ,Trunc(DBarrier), Dx, by eae 
Tr (Done) THEN 
Exec(’c:\command.com’,’/c copy done.bat answer.bat’); 
END. 
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B. SUBROUTINES 
1. SignalExcess 
iat Signald; 


{ Determines signal excess for the receiver and the 
target based on the bistatic active sonar equation. } 


Interface 
Weast CC = &.665889636; 


Pemer1on SignalExcess (Sx, Sy, Tx; Ty, Rx; Ry, PL, a, 
Tome Loman, Nii cie AG, DT =: Real): Real; 


Implementation 


Beeeeton otgnalexcess(Sx, Sy, Tx, Ty, Rx, Ry, PL, a, 
TSmax, TSmin, NL, SL, AG, DT : Real): Real; 


ware THR, ThS, Rst, Rtr, Rsr, an, TS : Real; 


Memeo n Lust, TLtr, Ftemp >: Real; 
BEGIN 
Rst := SORT(SOR(Tx - Sx) + SOR(Ty - Sy)); 
Rtr := SORT(SOR(Tx - Rx) + SOR(Ty - Ry)); 
Poon — OR (SOR (R= SX) + SOR(RY = Sy)); 
mies Tx — Sx; 
ne = Ix — Rx? 
Ite (kRst + RtEr) < (Rsr + 4800*PL/3000)) THEN 
SignalExcess := -100.00 
Bio 
BEGIN 


petargee Strength pare, TS} 


Vesey koe  Rerly *{nli*n2 + 0.25*(SOR(Rsr) - 
SOR(Rst - Rtr))); 
o2e- ee oma — TSmam) + (TSmain*an) ; 


MeiacacMitootom loss part for both paths, TLst, TLtr} 


Mice — me tN ORSt~1000)) +0(a*Rst ) ; 

ines — se “LNiktyr*= 1000)) + (a*ReEr); 

SignalExcess := SL-TLst-TLtr-NL-DT+TS+AG; 
END ; 


END; 
END. 


aS 


2. NormalValue 
Unit Normal: 
{ This function determines the cumulative normal 


distribution for a given value by using a polynomial 
approximation. } 


Interface 
Function NormalValue(x : Real): Real; 
Implementation 
Const dil = 4298673478072; 
d2 = 2.11410061E-02; 
Qs = 3.27 762658 -0-c- 
a4 = 3,80036E-05-; 
a5 = 4.88906E-05; 
G6 =" 5-26 5n—-0 0" 


Bee = 05000 0 i 
Function NormalValue(x : Real): Real; 


Var px eer cere ee eed. 


BEGIN 
ce = OORT Cookie. a 
Le (Xt bes) “THEN 
Be > “ee 10 5 
EiGSe 
BEGIN 
Xl := 1 + dl*¥xt + d2*SOR(xt) + d3“EXP(o i.e 
d4*SOR(SOR(xt)) + d5*EXP( 5" fue a eee 
d6*EXP(Gc*lmGeere 
x2 S00 257 EXP ae be at ee 
EN 
Le Ax = OO pe THe 
OP, aera a. <7 
PoE 
xe so = xe 
NormalValue := px; 
END: 
END. 
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cme FormatG 


Unit Gamsform; 


{ hits eprecealre converts the output from the main 
program into a form which the GAMS linear program solver can 
read in.} 


Interface 


Procedure FormatG(WBarrier, DBarrier, Dx, Dy, Dz 
iit eqeit): 


Implementation 


Procedure AdvanceData(Var Cl : Integer; Dz, C : Integer; 
Voit aOueOLD  : Text }’; 


Pome, ou : Integer; 
at >: Real; 


BEGIN 
Ul Os ie le ae 
Bekete;,= 1. to dU DO 
Read (DataOutOLD,T); 
Ge] i * Dz; 
END ; 


Procedure FormatG(WBarrier, DBarrier, Dx, Dy, Dz 


Integer) ; 
Yar DataOutOLD, DataOutNEW etext: 
eee ae Ls 8C2 C3  aeeqenx - 
mee CB, CC, 2 > dnclecer > 
Temp > istevel 
Ton wllge 2 Sie sia eye 
BEGIN 


Potomac aOueOLD, F< \THES OUT’) ; 

Assign (DataOutNEW, 'F:\THES.DAT’ ); 

Rewrite (Dat aOutNEW) ; 

Write (DataOUTNEW, ‘SET 7s oe 

Se] Oe 

ees = 0 
: OF: 

=| 1 

WHILE (C3 <= WBarrier) DO 

BEGIN 


Se 
es 
II 
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Le (Ce 3 — ees iy 
Write (DataoutNEw, 00 eer 


ELSE 
BEGIN 
1S eS ees Sales) 
fp @ ue mena 0) 
EESe 
BEGIN 
LE 4@3 = 100) rae 
pe =" ae 
BBSE 
Al>-4 ole a 
END; 
Write (DataOutNEW, °, ’7 Exe eee 
Ze a ee 
LP’ (Zee, 1s) Ae 
BEGIN 
Writeln(DataOutNEW) ; 
Write (DataOUTNEW, ' 
Za. oe 
END; 
END 
So =e ee 
END? 


Write(DataOutNEW, ’/'); 
Writeln(DataOutNEW) ; 


Le =e 
Write (DataOutNEW, ’ Youre 
WHILE (C2 <= DBarrier) DO 

BEGIN 


IF (C2 = 0) THEN 
Write (DataOutNEW, ‘'00‘,C2) 


ELS 
BEGIN 
Le fe? 2 vO ae 
ype = 00 
ELSE 
BEGIN 
Te (C2 = 21000 rae 
ye 3 = ae 
Boor 
Typ := °' 
EMD: 
Write (DataOutNEW,’,‘’,Typ,C2); 
Zp =e Zee 
LP (Z S-13)° AEN 
BEGIN 
Writeln(DataOutNEW) ; 
Write (DataOUTNEW, ’ 
i — ell, 2 
END; 
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END 
C7 e727 = by: 
END > 
Write (DataOutNEwW, ‘/’); 
Writeln(DataOutNEW) ; 


Th ele 
Write (DataOutNEw, ’ ay a ee 
WHILE (Cl <= WBarrier) DO 

BEGIN 


Ie (Ca = 10S Sab slay 
Write (DataOutNEW, C1) 


Jo ase a 
BEGIN 
(Wietee(DacaOueNeW, »°,Ci) ; 
fe = 0 tel 
END: 
Peet 2 aloe Teh 
BEGIN 


Writeln(DataOutNEW) ; 
Write (DataOUTNEW, ’ ay 
(Se — ae 
IB 
Giese = Cl + Dz: 
END; 

Write (DataOutNEW, ‘'/;'); 

Writeln(DataOutNEW) ; 

Writeln(DataOutNEW) ; 


Writeln(DataOutNEW, ' TABLE Pea sOmeioe,Y ,P) '); 
Writeln (DataQOutNEW) ; 
Wa ¢= Trunc (WBarrier/ DZ); 
ein: —= (- 
ee = |: 
mer <= 1; 
ae 5 i 
on 0: 
Write (DataOutNEw, ’ al 
REPEAT 
Cl — Oe 
G2 -= 0: 
Co ete). 
WHILE {{(8FTCA MOD 8) <> 0) AND (C <= LZ)) DO 
BEGIN 
Write (DataOutNEW, TC1:7) ; 
mOe = TCA + i; 
Go 3G =. i: 
AB@e =a Gi + Dz 
END: 


Writeln(DataOutNEW) ; 

Reset (DataOutOLD) ; 

WHILE (C3 <= WBarrier) DO 
BEGIN 


aS 


WHILE (C2 ==“BBarrie: jee 


END; 


Writeln 


BEGIN 
TP aes <- PO) ur eh 
2 gee — we (00k 
BoE 
BEGIN 
£E {(e3"—) 100 Sra 
xp =" 0M 
ELSE 
ee: cae 
END; 
Le Vee LO ere 
yo s=. 008 
E hier 
BEGIN 
TE (C2 = 200) TEN 
Eo ae) 
ELSE 
lo = 
END: 
Write (DataQutNEW, Txp, C3,” . 7 fee ie 


LE GT = 1) see 
AdvanceData(Cl, Dz, CT, DataOutOLD) ; 
IGS Sas — a ie 
WHILE ((Cl <= WBarrier) AND @I{TCBS MOD ci. 
O}) Be 
BEGIN 
Read (DataOutOLD, Temp) ; 
Write (DataOutNEW, Temp:7:3); 


Che = Cle ae 
ECBat: = Tebwee ai: 
END: 


Readln(DataOutOLD) ; 
Writeln(DataOutNEW) ; 
GIy a= 10 

GZ. s= "C2 ee by; 


END- 
=e Oh 


G3, GE: Dee 


DataOutNEw) ; 


Writeln(DataOutNEW) ; 


There <S 


iz 


Pe, (C222 [Ze eel 
Write (DataOutNEw, ’ ‘: Vs 


Cre = 


CT + 1+ 


UNRil tC =) Ze 
Write(DataOutNew, ’;’); 
Close (DataOutNEW) ; 


END 
SB 
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APPENDIX B 


A. GAMS FORMULATION PROGRAM 


eve ETE Bistatic Sonar Problem for Thesis 
oli TLE By LT Aasgelir Gangsaas USN 
*——------------ GAMS AND DOLLAR CONTROL OPTIONS-------------- 


SOrPPUPPER OFESYMLIST OFFSYMXREF 


SeETLIONS 
PIMcCOL=0 , LIMROW =0 Poorer ier) DECIMALSG=3 
Pao iM— VO TTERLIM=10000, OPTCR =) 0 Sia 6 =314aq 
05 SS SS [semen ean Mi Glee te a i ee 


POSITIVE VARIABLE 
PR(X,Y) Gikcmialioie LON fel eot bormt- 


VARIABLE 
V value of the game; 


FQUATIONS 
@O1L (P) 
PROB; 


* maximize 


Pesub ject to 
SOn(P).. SUM tee ery) @PAYORPE (xX, Y,P)) =G= V; 
PROB. . SUMS weit xk, Vv) > =B= 15 


MODEL ZEROSUM /ALL/; 
SOLVE ZGEROSUM USING LP MAXIMIZING V; 


COL.M(P)$(COL.M(P) EQ EPS) = 0.0000; 
MieorPhAY V.b, PR.L, COL.M; 

file outputl /D:\TP\WORK\GAMOUT.DAT/ ; 
ete Outcpuel : 

mae Vol 6:3 /; 
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LOOP CA oie, 

put x.TL:/, y.TL:7, PROG go ee 
ng 
file output2 /D:\TP\WORK\GAMDOUT.DAT/, ; 
DUE VOUeEpUEZ.. 

LO@GP« Ps 

DUE D.TL 77s COle Mr ae 

ie 


B. FINAL OUTPUT CONVERSION 


Program CumResults; 


{ This procedure takes the output from the GAMS program 
and outputs the results to files to cumulate results.} 


Uses DOS, Normal, Signal4; 


COnSteE ese si200 01s: 
Sigma = 8; 

Var UNDATAL OINDATAZ TNDATAS : Tests 
INDATA4, OUTDATA, OUTDATATEMP = Tex. 
OUTDATAI : exe: 
SOV, LA El CPZ Breadsd * Real 
SlmTPeMp,. <a; oh eye, BR : Real; 
Temp, SE, TSmax, TSmin, NL : Real; 
Si, 2AG,~ Dt : Real; 

Nec ees bee : Integer; 
PulseL, TSpeed, Freq : Integer; 
DBarrier, Pingin. Loeal Time : Integer; 


{RR RERRERKRE EK RRR ER ER RRR KER RRR RRR eee 


Procedure ReadIn; 


{ This procedure reads in from a file (Thes.in) which 
contains information about the receiver, target and source 
needed for the program. } 


Var Datailn =: Text - 


BEGIN 
Assign (Datain, ’D: \TP\WORK\ THES 2h) 
Reset (DataIn) ; 
Readln(DatalIn) ; 


Readln(DatalIn, PulselL) ; 
Readln(DatalIn) ; 
Readln(DataIn, Freq); 
Readln(DatalIn) ; 
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Readin(DatalIn, PingIn) ; 
Readin(Dataln) ; 
Readln(DatalIn) ; 
Readin(DatalIn) ; 
Readin(DataIn, DBarrier) ; 
Readln(Dataln) ; 
Readln(DatalIn); 
Readln(DatalIn) ; 
Readin(DataIn, TSpeed) ; 
Peadin(Dataln) ; 
Readlin(DatalIn) ; 
Readin(Dataln) ; 
Readln(DatalIn) ; 
Readln(DataIn) ; 
Readin(DatalIn, TSmax) ; 
Readln(Dataln); 
Readin(Datain, TSmin) ; 
Readin(DatalIn); 
Readin(DataIn, NL) ; 
Readln(DatalIn) ; 
Readin(DatalIn, SL); 
Readin(DataIn); 
Readin(DatalIn, AG) ; 
Readin(DatalIn) ; 
Readlin(DatalIn, DT); 
END; { Procedure ReadIn} 


ee ERE ER EER Malin Program ie A A EK EEE REE HE ® 


BEGIN 
ReadIin; 
Sumremp := 0.0; 
Assign(INDATAI1, ‘D:\TP\WORK\GAMOUT. DAT’ ) ; 
Assign(INDATA2, 'D:\TP\WORK\GAMDOUT.DAT’ ) ; 
Assign (INDATA3, ‘D:\TP\WORK\THESSRCE.IN’ ); 
Assign (INDATA4, ‘D:\TP\WORK\NUM. TXT’ ) ; 
Assign(OUTDATA, 'D: \TP\WORK\OUTPUT. DAT’); 
Assign (OUTDATAL, ’D: \TP\WORK\OUTPUT1.DAT’):; 
Assign (OUTDATATEMP, ‘D:\TP\WORK\TEMP.OUT’ ) ; 
Reset (INDATAL1) ; 
Reset (INDATAZ2 ) ; 
Reset (INDATA3) ; 
Reset (INDATA4) ; 
Rewrite (OUTDATATEMP) ; 


Bmeqgoqd := SOR(Freq/1000) ; 

eee Ome EB regqsa/ (1 + Freqsq)}) + (40*Freqsq/(4100 + 
FreqSq)) + ((2.75E-04) *FreqSq) ; 

motralTime >= Trunc( (DBarrier/ (TSpeed*2))*60); 


Readin(INDATAI1,I); 
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Readln(INDATA4 ,N) ; 
Readln(INDATAS3) ; 
FOR C 22) 1 Ge aN) eee 
Readln(INDATA3,No, Sx, Sy) ; 
he =e ee Tees 
BEGIN 
Rewrite(OUTDATA) ; 
Rewrite (OUTDATAL1 ) 
END 
Jaq So 
BEGIN 
Append (OUTDATA) ; 
Append (OUTDATAL) ; 


END; 
Writeln(OUTDATA, ’ Source Position (in thousands o£f 
Yards). 3s 
Writeln(OUTDATA, ’ X: ',Se:5:1,° "7 Yo eee 
Writeln(OUTDATA) ; 


( 
( 
Writeln(OUTDATA,’ Value of the Game: ‘,1:5:3); 
Writeln(OUTDATA) ; 

( 


Writeln(OUTDATA,'’ Optimal mixed strategy (for the 
Receiver) '); 
Writeln(OUTDATA, ’ Xx vy probability a 
WHILE not eof (INDATA1) DO 
BEGIN 


Read! nt INDAT AI Rx Ry blk); 
TF (Pi > EPs) THEN 
BEGIN 
Writeln(OUTDATA,Rx:-6:1,RY 628-132. 
Writeln(OUTDATATEMP, Rx, Ry, P1) 
Eo 
END; 
Close (OUTDATATEMP) ; 
Writeln(OUTDATA) ; 
Writeln(OUTDATA,’ Dual Solution (Optimal Strategy, Eomeea- 


Targee) ae 
Writeln(OUTDATA, ' X probability am 
WHILE not eof (INDATA2) DO 
BEGIN 


Readln(INDATA2,P, P2); 
Reset (OUTDATATEMP) ; 
|, ee — 1a) 
IF (P2 > EPS) THEN 
BEGIN 
Writeln(OUTDATA, P25, 22: Zao 
WHILE not eof (OUTDATATEMP) DO 


BEG 
a ge 
Thy = Beater ei. 
Readin (OUTDATA TEMES Rx, Ry, Ei 
5 le oO 
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Coen ie TOtal lame) | DO 
BEGIN 
Seepemc lOmcbeeeescerSs, SY, Rx, Ry, Tx, Ty, 
PulseL, a, TSmax, TSmin, 
NIG Sib Vee ean ee 
ee ee Come rr 
Temp := 1 - NormalValue(SE/Sigma) 
EOSE 
Temp := Temp*(1-NormalValue(SE/Sigma) ); 
ie = Tyee hopeea, (5°10) )*Pangin: 
5. ae DRE abate Jian 
END; 
Smet =—moummbems {i = Temp) *Pi*P2 : 
IND 
NB) 
JEAN ee 
Writeln(OUTDATA) ; 
Writeln(OUTDATA, ‘Counterdetection probability: ’ 
PSum nema: o 23) 
Writeln(OUTDATA) ; 
Writeln(OUTDATA, eee Dee en NN KKK KR AK ORR AE * | . 
Temeeim OUTDATAIT, (N = 1):4,1:7:3,SumTemp:7:3); 
Writeln(N - 1); 
Close (OUTDATAI1) ; 
Close(OUTDATA) ; 
END . 
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